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ABSTRACT 

The  combustion  or  afterburning  of  fuel-rich  rocket 
exhaust  with  the  atmosphere  may  result  in  large  infrared  radiation 
emissions  which  can  play  a  significant  role  in  the  design  of 
missile  base  components  and  missile  defense  systems.  Current 
engineering  level  models  neglect  turbulent-chemistry  interactions 
and  typically  underpredict  the  intensity  of  plume  afterburning 
and  afterburning  burnout.  To  evaluate  the  impact  of  turbulent- 
chemistry  interactions,  an  assumed  pdf  model  was  applied  to 
missile  plume  simulations  of  a  generic  booster.  Simulation 
results  reveal  turbulent-chemistry  interactions  to  have  a  large 
impact  on  plume  signatures  as  afterburning  burnout  was 
approached. 

I.  INTRODUCTION 

The  afterburning  of  missile  exhaust  with  the 
atmosphere  may  result  in  large  infrared  radiation  emission 
which  can  be  a  major  contributor  to  the  heat  transfer  to  the 
missile  base  (Reijasse  and  Delery,  1994).  As  altitude  is 
increased,  the  radiative  component  of  the  heat  transfer  rate 
will  eventually  show  a  large  drop  in  magnitude  (Kramer, 
1970)  resulting  from  the  shutdown  or  cessation  of 
afterburning  in  the  plume  (see  Fig.  1).  Shutdown  in  this 
context  does  not  refer  to  the  termination  of  the  missile 
engine,  but  to  the  cessation  of  the  combustion  taking  place 
between  the  missile  exhaust  and  the  atmosphere,  occurring 
with  continuous  engine  operation.  The  total  amount  of 
time  this  afterburning  occurs  and  its  rate  of  decay  during 
the  shutdown  event  will  determine  the  total  heat  transfer  to 
the  body  and  establish  a  design  criteria  for  the  components 
in  the  base  region.  The  radiative  emission  during  the 
shutdown  event  also  has  important  implications  for  the 
development  of  missile  defense  systems.  Consequently,  a 
need  exists  to  accurately  characterize  missile  plume 
afterburning  shutdown  or  cessation  events. 

The  character  of  afterburning  shutdown  events 
has  been  observed  to  vary  among  different  propulsion 
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exhibit  a  gradual  drop  in  total  radiant  intensity  over  a  wide 
altitude  range  while  others  shutdown  very  rapidly  over  a 
much  narrower  range.  Application  of  engineering  level 
modeling  techniques  (Dash,  et  ah,  1980)  to  predict 
afterburning  cessation  has  been  relatively  successful  for 
systems  that  exhibit  the  gradual  drop-off  type  of  shutdown 
event.  However,  for  systems  showing  abrupt  shutdown, 
engineering  models  have  been  less  successful. 

To  improve  understanding  into  the  physical 
mechanisms  producing  afterburning  shutdown,  Calhoon 
(2000)  investigated  the  shutdown  characteristics  of  a 
generic  amine  booster  system  within  the  framework  of  a 
computational  parametric  study.  Several  mechanisms  were 
investigated  to  explain  the  shutdown  behavior  of  this 
system  which  included:  1)  a  plume  shear  layer 
relaminarization  phenomenon,  2)  a  Damkohler  number 
effect  and  3)  a  classical  flame  extinction  mechanism.  The 
relaminarization  mechanism  was  found  to  be  implausible 
because  plume  shear  layer  Reynolds  numbers  were  well 
above  transition  after  shutdown  had  occurred.  The 
Damkohler  number  mechanism  was  found  to  be  plausible 
and  indeed  the  only  shutdown  mechanism  modeled  within 
most  commercially  available  codes.  This  results  from  the 
assumptions  used  to  model  the  mean  reaction  rate  in  the 
species  conservation  equations.  The  Damkohler  number 
mechanism  was  found  to  be  responsible  for  the  gradual 
drop-off  rates  of  plume  radiation  produced  by  most  codes. 
The  classical  flame  extinction  mechanism  was  found  to  be 
a  possible  explanation  for  the  observed  rapid  shutdown 
behavior  of  some  systems. 

The  flame  extinction  phenomenon  investigated  by 
Calhoon  (2000)  is  a  result  of  the  interaction  of  turbulence 
and  chemistry  at  small  spatial  scales.  The  model  used  by 
Calhoon  (2000)  to  account  for  this  phenomenon  was  an 
elementary  one  applicable  to  high  speed  compressible 
flows.  Other  models  do  exist  to  more  generally  account 
for  turbulent-chemistry  interactions  which  are  applicable  to 
compressible  flows.  Among  the  most  comprehensive 
techniques  are  the  assumed  probability-density-function 
(pdf)  method  (Gaffney,  et  al.,  1994)  and  the  compressible 
extension  of  the  pdf  evolution  equation  method  (Hsu  et  al.. 
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1994,  Delame  and  Pope,  1998).  Though  providing  a  more 
accurate  description  of  the  higher  order  statistics  of  the 
turbulent  scalar  fields,  the  pdf  evolution  equation  method 
has  not  been  shown  to  yield  significantly  better  results  than 
the  assumed  pdf  method  when  applied  to  compressible 
flows.  The  pdf  evolution  equation  method  is  also 
computationally  expensive  and  may  become  intractable 
when  applied  to  flows  including  shock  waves  (Delaine  and 
Pope,  1998).  The  assumed  pdf  method,  on  the  other  hand, 
is  computationally  inexpensive  and  offers  a  viable 
approach  to  accounting  for  turbulent-chemistry 
interactions  in  complex,  large  scale  flows  of  practical 
interest.  However,  the  assumed  pdf  method  does  not 
directly  account  for  strain  extinction  effects  and  may  not 
accurately  capture  this  phenomenon. 

For  high  speed  flows  of  interest  here,  turbulent- 
chemistry  interactions  have  been  shown  to  enhance 
burning  (Gaffney  et  al.,  1992)  for  flames  far  from 
equilibrium.  Including  turbulent-chemistry  interactions 
within  rocket  plume  simulations  has  the  potential  to 
improve  predictive  capability  because  engineering  level 
modeling  techniques  typically  underpredict  afterburning 
plume  temperatures  and  emissions  near  the  shutdown 
regime.  Including  turbulent  fluctuations  within  the 
radiation  calculations  may  also  significantly  impact 
predictions  (Pearce  and  Varma,  1981).  However,  this 
aspect  of  the  problem  was  no  considered  in  this  study.  The 
effect  of  turbulent  fluctuations  was  only  considered  for  the 
flowfield  predictions. 

The  objective  of  this  study  was  to  assess  the 
impact  of  turbulent-chemistry  interactions  on  the 
afterburning  and  afterburning  shutdown  characteristics  of  a 
generic  missile  system.  This  assessment  will  provide 
guidance  for  further  investigation  and  for  the  enhancement 
of  engineering  level  models.  This  work  is  a  continuation 
of  a  previous  study  (Calhoon,  2000)  which  investigated 
mechanisms  influencing  the  afterburning  shutdown 
characteristics  of  rocket  exhaust  plumes. 

In  the  following  sections  a  brief  review  of  the 
assumed  pdf  turbulent-chemistry  interaction  model  is  first 
given.  The  computational  methodology  used  within  the 
simulations  is  then  presented  followed  by  a  description  of 
the  generic  missile  geometry  and  engine  model  used. 
Results  are  then  presented  for  the  prediction  of 
afterburning  shutdown  for  this  missile  configuration  using 
the  assumed  pdf  model  followed  by  conclusions  that  may 
be  drawn  from  the  study.  These  result  indicate  that 
turbulent-chemistry  interaction  do  play  a  significant  role  in 
the  afterburning  and  afterburning  shutdown  characteristics 
of  missile  exhaust  plumes. 

II.  ASSUMED  PDF  METHOD 

The  effect  of  turbulent-chemistry  interactions 
within  the  context  of  Reynolds  averaged  Navier-Stokes 
simulations  (RANS)  appears  as  additional  unclosed  terms 
in  the  governing  steady  state  conservation  equations. 


These  unclosed  terms  include:  1)  the  time  averaged 
reaction  rate  term  wk  in  the  species  conservation 
equations  and  2)  temperature-species  correlations  in  the 
state  equations.  The  reaction  rate  is  a  highly  nonlinear 
function  of  temperature  and  species  concentrations  and  its 
time  average  may  be  expressed  as, 

_  oooo  1  1  1 

JJJ-I'J  wk  (A  T.  Yk  )P(p,  T,  Yk  )dYl...dYKdTdp  (1) 

0  0  0  0  0 

where  T,  p,  Yk  and  wk  are  the  temperature,  density  and  /fth 
species  mass  fraction  and  reaction  rate  for  species 
l<k<K.  P(p,T,Yk)  is  the  single-point  joint  probability 


density  function  of  p.  T  and  Yk  and  represents  the  combined 
effects  of  turbulent  transport,  both  large  and  small  scale, 
and  molecular  diffusion.  The  time  averaged  state  equation 
for  a  mixture  of  perfect  gasses  may  be  written  as, 


P  =  pR°X 


L  <T><  Y,  > 


k-\  Wk  k= 1  rr  k 

where  the  brackets  (<  >)  and  tilde  represent  conventional 
and  Favre  time  averaging,  p  is  the  pressure,  R°  is  the 
universal  gas  constant  and  Wk-  is  the  Mi  species  molecular 
weight.  The  fluctuating  components  T"  and  Y"  of  T  and 


+  R°P  I 


£ <T"Y”> 


Wt 


(2) 


Yk,  are  defined  by  T  =<T  >  +T"  and  Yk  =<  Yk  >  +Y" , 
respectively. 

Temperature-species  correlations  similar  to  what 
is  seen  in  Equ.  (2)  also  appear  in  the  time  averaged  caloric 
equation  of  state,  <  h  >=<  e  >  +p  /  p ,  relating  the  internal 


energy  to  enthalpy.  Expressing  the  species  specific 
enthalpies  in  terms  of  the  standard  polynomial  curve  fits, 
the  time  averaged  enthalpy  becomes 

<h>=  f(pT"Y",...,pT"NYk)  where  N  is  the  number  of 
coefficients  in  the  fits. 

To  close  the  governing  flow  equations,  an 
approximation  for  the  pdf  P  must  be  specified.  Within  the 
context  of  the  assumed  pdf  approach,  the  shape  of  P  is 
specified  in  terms  of  functions  which  may  be 
parameterized  by  their  lower  moments.  The  first 
approximation  used  within  this  approach  is  to  assume 
statistical  independence  of  p,  T  and  Yk  so  that  P  may  be 
expressed  as, 

P(p,TJk)=Pp(p)PT(T)PQ(Yk)  (3) 


where  P0  PT  and  Pq  are  the  marginal  pdf  s  of  p,T  and  Yk, 
respectively.  This  assumption  is  questionable.  However, 
assumed  pdf  methods,  cast  in  terms  of  a  mixture  fraction 
and  chemical  progress  variable  formulation,  have  been 
shown  to  produce  good  results  for  a  range  of  low  speed 
turbulent  reacting  flow  problems  using  this  assumption. 
Consequently,  this  assumption  is  carried  over  to  the  more 
general  formulation  described  here.  This  assumption  also 
results  in  <T"Y">=  0  since  the  variables  are  not 
correlated. 

Baurle  and  Girimaji  (1999)  investigated  a 
modification  to  this  procedure  which  relaxed  this 
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assumption  and  showed  potential  for  improvement. 
However,  further  research  is  required  to  resolve 
realizability  issues  that  were  identified.  Therefore,  this 
modified  approach  was  not  considered  further  in  this 
study. 

Consistent  with  work  of  Gaffney,  et  al.  (1994), 
the  marginal  pdf  s  P/h  PT  and  Pq  are  specified  as  follows. 
A  delta  function  is  assumed  for  the  marginal  pdf  of  density 
so  that  Pp=b(p  -  p).  For  the  temperature  pdf PT,  a  beta 


function  was  used.  The  beta  function  has  been  shown  to 
accurately  capture  scalar  mixing  in  homogeneous 
turbulence  (Girimaji,  1991).  FV  was  taken  as, 


with. 


and, 


9p'~l  (1-0)  A-1 

prV)=  '  HA+A) 


P,=<6> 


P2  =(!"<*>) 


T-T  . 

q  _  min 


r(AW2) 

<  9  >  (1-  <  9  >) 


-1 


<  9"9"  > 

<  o  >  (l-  <  e  >) 
<e"0"> 


-l 


T  —T 


<0>=<T>  Tmm 


T  -T 

max  mm 


<  0  "0 "  >=  - 


<rr> 

(T  -T  ) 

V  max  mm  / 


(4) 


(5) 


(6) 


where  F (x)  is  the  gamma  function.  The  variance  <  9" 9"  > 
was  explicitly  limited  by  the  realizability  constraint, 

<9"9"  ><<9>(\-<9>)  (7) 

The  temperature  limits  Tmin  and  Tmax  were  specified  as, 

Tmin  =  max|T0 ,  T  -  n^<T"2  >]  ^ 

^inax  =  min[7’1,7’  +  iv\l<  T"2  >] 


where  the  coefficient  n  was  specified  to  have  a  value  of  3 
following  Gaffney,  et  al.  (1992).  T0  and  77  are  problem 
dependent  absolute  minimum  and  maximum  realizable 
temperature  limits.  Specifying  Tmin  and  Tmax  in  this  manner 
gives  the  pdf  a  symmetric  distribution  for  small 
temperature  variance  while  allowing  for  nonsymmetric 
shapes  when  the  limits  are  capped  by  T0  and/or  7). 
Allowing  nonsymmetric  pdf  shapes  also  allows  for  the 
physically  realistic  situation  of  temperature  intensities 


being  greater  than  1  (i.e.,  v<  T"2  >  IT  >  1 ). 

For  the  pdf  Pq,  a  multivariate  beta  distribution 
was  used  and  has  the  following  form, 

p  (Y  Y  Y  x- F(A +A +-  +  ^) 

rQ\J  15  1  2  ’***’  1  K  ) 

""  '  (9) 


r(A)r  (P2)...r(pK) 


with. 


Pk  =<Yk> 


P-S 

~ Q~ 


(10) 


where, 

S=fj<Yk>2,Q=fj<Y£Y£>  (11) 

k—\  k-l 

The  quantity  Q  is  the  sum  of  the  species  variances  and  is 
called  the  turbulent  scalar  energy.  As  with  PT  ,  a 
realizability  condition  must  be  explicitly  imposed.  This 
condition  was  Q  <  (1  -  S) . 

With  Equs.  (4)  -  (11)  the  joint  pdf  of  p,  T  and  Yk 
is  completely  specified  given  <T>,  <T"T">,  <Yk  > 

and  Q.  The  mean  reaction  rate  may  then  be  calculated 
from  Equ.  (1).  Numerically  integrating  Equ.  (1)  can 
become  very  computationally  expensive  for  a  large  number 
of  species.  Fortunately,  due  to  the  form  of  pdf  s,  integrals 
involving  species  and  Po  can  be  evaluated  analytically 
(Gaffney  et  al.,  1994).  The  integrals  involving  T  and  PT 
were  evaluated  numerically  and  stored  in  look-up  table 
form  as  a  function  of  T  and  the  variance  of  T  ( Var(T )) 
(e.g.,  Gerlinger  et  al.,  1999). 

The  remaining  issue  with  respect  to  the  evaluation 
of  Equ.  (1)  is  the  specification  of  the  temperature  variance 
and  scalar  energy.  These  quantities  were  evaluated  using 
modeled  evolution  equations  derived  from  the  species  and 
energy  equations  and  are  given  by  (Gaffney,  et  al.  1992; 
Baurle,  et  al.,  1994a), 


Dp<g> 
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(13) 


K . 


+2Z« 

k=\ 


where  e  is  the  internal  energy,  <g>  is  the  internal  energy 
variance,  Uj  is  the  velocity,  y  is  the  ratio  of  specific  heats 
and  k  and  s  are  the  turbulent  kinetic  energy  and  energy 
dissipation  rate.  The  coefficients  v,  vT,  Sc  and  Pr  are  the 
kinematic  viscosity,  turbulent  eddy  viscosity,  and  Schmidt 
and  Prandtl  numbers,  respectively.  The  length  scales  of 
the  temperature  and  species  fluctuations  were  assumed  to 
be  proportional  to  the  velocity  length  scale  and  specified  in 
terms  of  the  turbulent  Schmidt  and  Prandtl  numbers,  ScT 
and  PrT  ,  respectively.  The  last  term  in  Equ.  (12)  is  a 
dilatation  term  resulting  from  compressibility  effects 


3 

American  Institute  of  Aeronautics  and  Astronautics 


AIAA-2000-3388 


(Gaffney,  et  al.  1992).  The  last  term  in  Equ.  (13)  is  a 
chemical  source  term  which  may  be  evaluated  using  the 
assume  pdf  s  as  described  by  Gaffney  et  al.  (1994). 

To  obtain  the  temperature  and  the  temperature 
variance  from  the  mean  internal  energy,  <e>,  and  its 
variance,  <g>,  temperature  and  species  fluctuations  were 
neglected  so  that  the  temperature  may  be  evaluated  in  the 
usual  manner  by  solving  <  e  >  ^  ><  et  >=  0  .  To 

obtain  Var(T),  the  specific  heat,  c,,,  was  linearized  about 
the  mean  temperature  so  that, 

<T”2  >«  (14) 

K 

where  cv  =  cv  (<  T  >,  <  Yk  >) . 

The  dissipation  coefficients  Cg  and  C()  in  Equs. 
(12)  and  (13),  respectively,  were  specified  as  Cg  =  Cq  = 
.546  by  matching  the  modeled  scalar  decay  rate  to  that 
obtained  from  the  DNS  results  of  Eswaran  and  Pope 
(1988)  and  McMurtry,  et  al.  (1993)  for  forced  isotropic 
turbulence.  The  production  coefficients  in  these  equations 
were  specified  as  Ct  =  .2  and  C2  =  .5  as  will  be  discussed 
in  a  later  section. 

III.  COMPUTATIONAL  METHODOLOGY 

The  missile  plume  flowfields  in  this  study  were 
analyzed  using  CRAFT  (Sinha,  et  al.,  1992)  which  is  a 
structured,  finite-volume  code  that  solves  the  compressible 
Navier-Stokes  (NS)  equations.  The  solver  is  fully  implicit 
and  uses  Roe/TVD  numerics  for  the  inviscid  fluxes  and 
second-order  central  differencing  for  viscous  and  diffusive 
terms.  CRAFT  is  sufficiently  general  to  simulate  finite- 
rate  chemistry  and  multi-phase  flows  and  includes  standard 
polynomial  curve  fits  for  the  thermodynamic  properties. 
The  code  includes  a  variety  of  advanced  turbulence  model 
formulations  (Kenzakowski,  et  al.,  2000).  However,  in 
this  study,  only  the  standard  k-s  model  with 
compressibility  correction  was  used.  CRAFT  also  includes 
an  implementation  of  the  assumed  pdf  turbulent-chemistry 
interaction  model  described  in  the  previous  section.  The 
CRAFT  code  achieves  good  computational  efficiency 
through  a  fully  parallelized  implementation  of  the  flow 
solver  using  a  combined  shared  and  distributed  memory 
(using  MPI)  data  parallel  model.  Computational  efficiency 
is  also  enhanced  using  a  parabolized  form  of  the  NS  flow 
solver  (PNS)  for  parabolic/hyperbolic  regions  of  the  flow. 
This  PNS  solver  includes  a  sub-stepping  option  to  allow 
for  axial  grid  refinement  in  regions  where  flow  properties 
change  rapidly. 

Radiation  calculations  were  decoupled  from  the 
flowfield  simulations  and  were  carried  out  as  described  by 
Ludwig  et  al.  (1981)  and  Nelson  (1987).  Within  this 
approach,  radiation  transport  equations  were  solved  using 
a  band  model  for  the  gas  phase  absorption/emission.  The 
calculations  for  the  present  study  used  a  wide  band  pass  to 
encompass  emissions  from  OH,  CO,  C02  and  H:0.  Total 


radiation  intensity  predictions  were  also  made  using  a  field 
of  view  large  enough  to  encompass  the  entire  plume. 

IV.  ROCKET  PLUME  SIMULATIONS 

The  simulations  considered  in  this  study  were 
carried  out  for  a  generic  axisymmetric  amine  booster 
described  in  detail  by  Calhoon  (2000).  The  exhaust  nozzle 
was  assumed  to  have  an  area  ratio  of  10  and  protrude  aft  of 
the  booster  base.  The  entire  missile  body-base -plume  was 
simulated  using  CRAFT’s  NS  and  PNS  solvers.  The 
elliptic  regions  along  the  body  and  in  the  base  region  were 
computed  using  the  NS  solver.  The  hyperbolic  plume  flow 
was  calculated  using  the  PNS  option.  Each  axisymmetric 
solution  used  49,000  grid  points  to  resolve  the  body-base 
region  and  700,000  points  for  the  plume  region.  Each 
grid,  both  body-base  and  plume,  was  manually  adapted 
around  regions  with  high  gradients  to  ensure  proper 
resolution  of  the  flow  features. 

The  inflow  boundary  conditions  at  the  engine 
nozzle  exit  plane  were  specified  from  a  separate  nozzle 
flow  calculation  described  in  detail  by  Calhoon  (2000).  A 
standard  9  species  -  10  reaction  step  mechanism  (Calhoon, 
2000)  for  H2/CO  oxidation  was  used  for  this  calculation  as 
well  as  for  the  plume  simulations.  For  this  nozzle 
simulation  the  flow  was  fairly  uniform  across  the  exit, 
except  in  the  boundary  layer  region  along  the  wall,  with  a 
core  temperature  of  ~  1330  K.  The  fuel-rich  conditions  of 
the  engine  were  evident  from  an  excess  of  CO  and  some 
H2  in  the  core  with  0_  virtually  depleted.  This  nozzle 
solution  was  assumed  to  be  invariant  with  respect  to 
altitude  and  was  used  for  all  the  body-base-plume 
calculations.  This  solution,  however,  did  not  account  for 
turbulent-chemistry  interactions  within  the  nozzle.  At  this 
inflow  the  temperature  intensity  was  specified  as  5%  while 
the  turbulent  scalar  energy  was  specified  as  0%. 

V.  RESULTS  AND  DISCUSSION 

The  body-base  and  plume  flowfields  were 
calculated  at  three  altitudes  (25  km,  30  km  and  35  km)  both 
with  and  without  the  pdf  model  described  earlier.  These 
altitudes  were  selected  because  they  span  the  afterburning 
shutdown  regime  for  this  generic  booster,  given  an 
assumed  trajectory  profile  (Calhoon,  2000).  The  Reynolds 
and  Mach  numbers  based  on  freestream  conditions  and  the 
body  radius  were  Rex  =  4x1 06,  2x1 06  and  lxl  06  and  Mx  = 
2.6,  3.2  and  3.9  for  the  altitudes  of  25,  30  and  35  km, 
respectively. 

Plume  simulations  were  carried  out  at  these  three 
altitudes  for  three  different  levels  of  modeling.  The  first 
assumed  turbulent-chemistry  interactions  to  be  negligible. 
This  approach  was  termed  the  ‘laminar’  reaction  rate 
model.  The  second  model  included  turbulent-chemistry 
interactions  via  the  temperature  pdf  only,  while  the  third 
approach  included  both  the  temperature  and  species  pdf  s. 
These  three  levels  of  modeling  were  considered  not  only  to 
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assess  the  effect  of  turbulent-chemistry  interactions  on 
missile  plume  afterburning  shutdown,  but  also  to  assess  the 
effectiveness  of  the  different  aspects  of  the  assumed  pdf 
model  itself.  For  the  pdf  cases,  the  absolute  temperature 
limits,  T„  and  T ,,  required  to  evaluate  the  pdf  temperature 
limits,  Tmin  and  Tmax,  in  Equ.  (8)  were  specified  as  follows. 
With  the  freestream  temperatures  for  the  specified  flight 
conditions  ranging  from  222  to  239  K,  the  absolute 
minimum  temperature  was  specified  as  T0  =  210  K.  The 
maximum  absolute  temperature,  Th  was  specified  as  2800 
K  which  corresponded  to  the  maximum  temperature  in  the 
boundary  layer  along  the  nozzle  wall  (Calhoon,  2000). 
With  these  absolute  temperatures,  the  maximum  realizable 

temperature  intensity  was  V<  T"2  >  IT  ~1.6,  which 
corresponds  to  the  condition  when  freestream  and  nozzle 
boundary  layer  fluid  co-exist  in  an  unmixed  state. 

As  described  by  Calhoon  (2000),  the  altitude 
conditions  of  25,  30  and  35  km  characterize  the 
afterburning  shutdown  regime  of  the  present  generic 
booster  system  when  simulated  using  the  ‘laminar’ 
reaction  rate  model.  For  this  model,  the  plume  at  25  km  is 
burning  vigorously,  while  at  30  km  it  has  progressed  deep 
into  the  shutdown  regime.  At  35  km  the  plume  flame  is 
almost  completely  extinguished.  Application  of  the 
assumed  pdf  model  to  these  flight  conditions  was  found  to 
have  a  substantial  impact  on  the  afterburning 
characteristics  of  this  system.  For  example,  Fig.  2  presents 
contour  plots  of  temperature  and  C02  mole  fraction  at  30 
km  for  both  the  laminar  rate  and  temperature  pdf  models. 
Simulations  using  both  the  temperature  and  species  pdf 
will  be  discussed  later.  From  this  figure  the  plume  is  seen 
to  rapidly  expand  due  to  the  highly  underexpanded 
condition  at  the  nozzle  exit.  The  plume  barrel  shock  is 
clearly  evident  in  the  temperature  field  along  with  its 
subsequent  reflection  off  the  axis  of  symmetry  and 
interaction  with  the  plume  shear  layer.  Further 
downstream  a  series  of  weaker  reflections  persist  until  the 
plume  shear  layer  merges  with  the  axis.  For  the  laminar 
rate  case,  plume  ignition  is  delayed  downstream  until 
approximately  three  times  the  distance  between  the  exit 
plane  and  the  barrel  shock  reflection  point.  This  is  evident 
from  the  figure  by  the  delay  in  the  rise  of  C02  in  the  plume 
shear  layer.  This  ignition  delay  also  results  in  a  long  delay 
in  the  plume  shear  layer  temperature  rise.  For  the 
temperature  pdf  case,  however,  it  is  clear  the  model  has  a 
significant  impact  on  the  combustion  processes  occurring 
within  the  plume.  For  the  pdf  case,  ignition  occurs  close  to 
the  missile  base  as  evidenced  by  the  rapid  rise  in  shear 
layer  C02.  Further  downstream,  the  plume  vigorously 
afterburns  resulting  in  higher  plume  temperatures. 

The  enhanced  burning  realized  for  the 
temperature  pdf  model  at  this  altitude  was  found  for  the 
other  altitudes  as  well.  Fig.  3  presents  a  comparison  of  the 
laminar  and  temperature  pdf  rate  model  predicted  mean 
temperatures  along  the  plume  axis  for  the  three  altitudes. 
In  this  figure  the  axial  coordinate  has  been  scaled  by  xrefj. 


This  length  scale  is  the  distance  between  the  nozzle  exit 
plane  and  the  barrel  shock  reflection  point  at  the  axis  of 
symmetry.  This  reference  length  was  found  to 
approximately  scale  the  results  with  respect  to  altitude  so 
that  a  fixed  value  of  x/xreli  corresponds  to  the  same  relative 
location  in  the  plume  for  each  case  (Calhoon,  2000).  From 
Fig.  3,  the  pdf  model  is  seen  to  enhanced  afterburning  at 
all  altitudes  as  evidenced  from  the  higher  peak 
temperatures  in  the  plume  farfield.  At  25  km,  when  the 
laminar  model  is  vigorously  afterburning,  the  difference 
between  the  two  models  is  relatively  small.  However,  as 
altitude  is  increased  and  afterburning  shuts  down  for  the 
laminar  model,  the  differences  become  more  pronounced. 
At  35  km  the  pdf  case  is  still  afterburning  while  the 
laminar  rate  case  is  nearly  extinguished  resulting  in  a 
temperature  difference  of  approximately  150  K  in  the 
farfield  plume.  As  altitude  is  increased,  the  pdf  model 
shows  a  gradual  migration  of  the  ignition  point  further 
downstream  consistent  with  the  laminar  rate  model. 
However,  the  temperature  pdf  model  has  delayed  the  onset 
of  afterburning  shutdown  significantly. 

Gaffney,  et  al.  (1992)  found  temperature 
fluctuations  to  enhance  chemical  reactions  and  ignition  for 
high  speed  H2-air  shear  layer  flames  that  were  far  from 
equilibrium.  The  same  observation  may  be  made  here  with 
regard  to  CO/H2-air  combustion  occurring  in  these  plume 
flows.  For  the  present  plume  application,  the  laminar  rate 
model  produces  plume  flames  that  are  far  from  equilibrium 
and  on  the  edge  of  burnout.  When  the  pdf  model  is 
applied,  temperature  fluctuations  enhance  product 
formation  resulting  in  higher  plume  temperatures  and  a 
delay  in  afterburning  shutdown.  This  may  be  seen  in  Figs. 
4  and  5  which  are  plots  of  C02  mole  fraction  and 
temperature,  respectively,  at  x/xrefi  =  3  for  each  altitude. 
The  transverse  coordinate  in  these  figures  has  been  scaled 
by  8  which  is  the  location  at  which  the  N2  mole  fraction  is 
99%  of  its  freestream  value  for  the  laminar  rate  cases. 
From  Fig.  4  the  temperature  pdf  model  is  seen  to  enhance 
the  formation  of  C02  for  all  altitudes  except  25  km.  At  this 
altitude  the  laminar  rate  model  produces  more  C02  than 
the  pdf  model.  However,  the  H20  formation  (not  shown) 
at  25  km  is  much  greater  for  the  pdf  model,  offsetting  the 
lower  C02  formation  so  that  the  pdf  model  produces 
higher  peak  temperatures  for  all  altitudes  (Fig.  5).  The 
enhancement  of  the  H20  chain  within  the  mechanism 
results  in  a  large  enhancement  of  OH  for  the  pdf  as  well 
(Fig.  6).  For  the  laminar  rate  cases,  Fig.  4  shows  a  large 
drop  in  peak  C02  at  this  streamwise  location  as 
afterburning  shuts  down  and  the  plume  ignition  point 
moves  further  downstream.  The  pdf  model  shows  a 
similar  trend  but  to  a  lesser  extent,  indicating  the  plume  is 
shutting  down  slower.  Also,  similar  to  Fig.  3,  Fig.  5  shows 
the  impact  of  the  pdf  model  on  the  temperature  field  to 
increase  with  increasing  altitude. 

The  large  impact  of  the  assumed  pdf  method  on 
plume  combustion  processes  results  from  the  influence  of 
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temperature  fluctuations  on  the  evaluation  of  the  mean 
chemical  reaction  rate.  Fig.  7  presents  a  plot  of  the 
temperature  fluctuation  intensity  as  a  function  of  altitude  at 
x/xreft  =  3.  The  fluctuation  levels  are  seen  to  be  quite  large 
at  this  streamwise  station  and  increase  in  magnitude  with 
increasing  altitude.  These  large  values  are  a  result  of  the 
excitation  of  the  plume  shear  layer  at  the  point  where  the 
reflected  barrel  shock  intersects  the  layer.  The  increasing 
magnitude  of  the  fluctuations  indicates  the  pdf  model  is 
becoming  increasingly  important  with  increasing  altitude  . 
Consequently,  neglecting  these  fluctuations  would  be 
inappropriate  especially  as  afterburning  shutdown  is 
approached.  The  primary  effect  of  the  temperature  pdf 
model  is  on  the  ignition  process  so  that  prior  to  shutdown 
the  model  is  much  less  significant  as  seen  for  the  25  km 
case. 

From  Fig.  7  the  model  predicts  the  temperature 
intensity  to  increase  with  increasing  altitude.  From  Equ. 
(12),  it  is  unclear  what  the  source  of  this  trend  is.  Fig.  8 
presents  a  plot  of  the  production,  dissipation  and  dilatation 
terms  of  Equ.  (12)  at  x/xreji  =  3  for  each  altitude.  From  this 
figure  the  magnitude  of  each  term  is  found  to  be 
decreasing  with  increasing  altitude.  Also  note  that  the 
dilatation  term  can  have  both  positive  and  negative  values 
and  is  small  compared  to  the  production  and  dissipation 
terms.  Consequently,  compressibility  effects  are  not 
explicitly  playing  a  significant  role  in  the  evolution  of  the 
energy  variance  in  Equ.  (12).  Flowever,  compressibility 
effects  are  implicitly  influencing  Equ.  (12)  through  k  and  s 
which  evolve  from  the  compressibility  corrected  k-s  model 
used  in  the  simulations.  Fig.  9  presents  the  sum  of  the 
production,  dissipation  and  dilatation  terms  in  Equ.  (12)  at 
this  same  streamwise  location.  From  this  figure,  the  total 
source  term  contribution  to  Equ.  (12)  is  found  to  be  rapidly 
decreasing  with  increasing  altitude.  This  indicates  that  the 
increasing  temperature  intensity  fluctuations  seen  in  Fig.  7 
are  a  result  of  upstream  influences  and  not  a  result  of 
enhanced  relative  production  within  the  shear  layer. 
Examination  of  contour  plots  of  turbulent  kinetic  energy 
and  temperature  intensity  show  the  plume  shear  layer  to  be 
greatly  influenced  by  the  reflected  barrel  shock.  This 
reflected  shock  (see  Fig.  2)  excites  the  shear  layer  causing 
the  layer  to  thicken  and  it  enhances  the  production  of 
turbulence  and  temperature  fluctuations.  This  trend  is 
consistent  with  the  work  of  Norris  and  Edwards  (1997), 
who  applied  a  much  more  general  unsteady  large-eddy 
simulation  technique  to  compute  high-speed  reacting 
exhaust  flows.  This  shock-shear  layer  interaction 
determines  the  turbulence  and  temperature  fluctuation 
levels  which  feed  into  the  evolving  shear  layer.  This 
interaction  is  the  apparent  source  of  the  trend  observed  in 
Fig.  7. 

The  large  increase  in  plume  temperatures  for  the 
temperature  pdf  simulations  correspondingly  results  in  a 
large  increase  in  plume  radiative  emissions.  Fig.  10  and  11 
present  comparisons  of  station  radiation  and  total  radiant 


intensity,  respectively.  From  Fig.  10,  the  station  radiation 
for  the  pdf  model  is  seen  to  be  significantly  larger  than  for 
the  laminar  rate  model.  Burning  for  the  pdf  model  at  35 
km  is  also  evident  while  the  laminar  model  is  nearly 
extinguished  at  this  altitude.  The  higher  station  radiation 
predictions  in  Fig.  10  produce  significant  differences  in  the 
total  plume  intensity  as  seen  in  Fig.  11.  As  noted  earlier, 
the  pdf  model  has  a  greater  influence  as  altitude  is 
increased.  ThistrendisalsoveryapparentinFig.il.  This 
figure  also  shows  the  pdf  model  to  delay  the  onset  of 
afterburning  shutdown  and  cause  the  plume  to  burn  to 
higher  altitudes  than  for  the  laminar  rate  model.  The  pdf 
model  also  significantly  changes  the  afterburning 
shutdown  rate,  slowing  it  over  what  is  seen  for  the  laminar 
rate  model. 

Though  changing  the  magnitude  of  the  emissions, 
the  pdf  model  did  not  change  the  character  of  the 
shutdown  event  seen  in  Fig.  11.  Both  the  laminar  and  pdf 
models  show  a  gradual  type  shutdown  behavior  which  is 
characteristic  of  a  Damkohler  number  effect  as  described 
by  Calhoon  (2000).  As  demonstrated  by  Calhoon  (2000), 
a  strain  rate  induced  extinction  phenomenon  may  be  the 
source  of  rapid  afterburning  shutdown  behavior  which  has 
been  observed  for  some  missile  systems.  A  strain  rate 
extinction  mechanism  is  a  consequence  of  turbulent- 
chemistry  interactions.  Flowever,  an  extinction  mechanism 
is  not  modeled  within  the  present  assumed  pdf 
formulation.  Therefore,  the  assumed  pdf  model  cannot 
capture  such  a  phenomenon. 

As  mentioned  in  Sec.  II.,  the  coefficients  C/  and 
C2  on  the  production  terms  in  Equs.  (12)  and  (13)  were 
specified  as  C;  =  .2  and  C2  =  .5.  The  value  C2  was 
specified  by  optimizing  the  coefficient  to  match  the  scalar 
fluctuation  data  of  Lockwood  and  Moneib  (1980)  for  a  low 
speed  nonreacting  jet.  Applying  the  same  optimization 
procedure  for  Equ.  (12)  resulted  in  a  value  of  .343  for  C’/. 
Figs.  12  and  13  present  the  sensitivity  of  the  temperature 
and  temperature  intensity  with  respect  to  variations  in  Cj  at 
x/x^  =  3  for  the  30  km  case.  From  Fig.  13  the  temperature 
intensity  is  seen  to  vary  by  ~  50  %  over  the  range  .1  <Cj< 
.343  while  the  peak  temperature  in  Fig.  12  varies  by  only  ~ 
6  %.  For  simulations  with  C/  =  .343.  the  temperature 
intensities  were  high  enough  to  force  the  plume 
afterburning  to  be  at  a  near  equilibrium  condition  for  each 
altitude.  Under  this  condition  the  normalized  C02  and 
HiO  mole  fractions  were  invariant  with  altitude  and 
afterburning  shutdown  was  not  initiated  within  the  25  to  35 
km  window  of  the  simulations.  This  result  was  unrealistic 
for  the  generic  booster  under  consideration  .  This  also 
indicates  that  the  low  speed  calibration  is  inappropriate  for 
application  to  the  high  speed,  highly  compressible  flows  of 
the  present  context.  The  present  value  of  C,  =  .2  was 
selected  because  it  was  close  to  the  calibrated  value  and 
allowed  afterburning  shutdown  to  be  initiated  within  the 
25  to  35  bn  window,  as  expected.  The  lower  value  of  Cj 
=  .1  showed  the  same  trends  for  mean  flow  quantities  as 
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for  the  present  value  of  Cy  =  .2,  but  to  a  slightly  lesser 
extent.  The  value  selected  here  for  C,  was  deemed 
appropriate  since  the  objective  of  this  study  was  to  make 
an  assessment  of  the  potential  impact  of  turbulent- 
chemistry  interactions  on  exhaust  plume  simulations,  and 
not  to  make  quantitative  predictions.  To  resolve  this  issue 
will  require  a  rigorous  validation  of  Equ.  (12)  for  the  high 
speed  flow  regime,  a  regime  in  which  little  highly  resolved 
experimental  data  exists  for  scalar  fluctuations.  In 
addition,  a  more  advanced  formulation  may  be  required 
which  includes  a  length  scale  equation  for  temperature 
fluctuations  (Chidambaram,  etal.,  1999). 

As  a  final  topic,  Fig.  14  presents  a  comparison  of 
the  mean  temperature  at  x/xrejj  =  3  calculated  using  the 
temperature  pdf  only  model  and  the  full  temperature  and 
species  pdf  formulation  described  in  Sec.  II.  From  the 
figure  it  is  apparent  that  the  full  temperature-species  pdf 
formulation  produced  virtually  identical  results  to  the 
temperature  pdf  only  model.  This  was  found  to  be  true  for 
the  other  flow  variables  as  well.  The  ineffectiveness  of  the 
assumed  species  pdf  model  to  significantly  change  the 
results  can  be  traced  to  the  prediction  of  the  scalar  energy, 
Q ,  from  Equ.  (13).  Fig.  15  presents  a  plot  of  Q  across  the 
plume  shear  layer  at  x/xreJi  =  3  calculated  using  Equ.  (13) 
both  with  and  without  including  the  chemical  source  term 
at  30  km.  Including  the  chemical  source  term  causes  Q  to 
be  driven  to  small  values  effectively  turning  off  the  species 
pdf  model,  producing  the  same  result  as  for  the 
temperature  pdf  only  model.  Eliminating  this  source  term 
from  Equ.  (13)  produces  large  values  of  scalar  energy  as 
seen  in  the  figure.  This  chemical  source  term  was  found  to 
be  dissipative  and  to  destroy  the  scalar  energy  and  drive  it 
to  small  values.  Baurle,  et  al.  (1994b)  observed  the  same 
behavior  for  a  high  speed  H2-air  flame.  In  another 
publication,  Baurle,  et  al.  (1995)  also  compared 
predictions  of  this  source  term  using  the  assumed  pdf 
model  and  a  more  comprehensive  pdf  evolution  equation 
model  for  the  same  jet  flame.  This  comparison  showed  the 
predictions  of  this  term  generally  to  be  of  opposite  sign 
and  to  be  different  in  magnitude  by  a  factor  of  five.  The 
present  results  support  the  assertion  of  Baurle,  et  al.  (1995) 
that  the  form  of  the  species  pdf,  Pq,  is  incapable  of 
computing  the  higher  order  moments  required  to  evaluate 
the  Q  chemical  source  term  with  any  reasonable  degree  of 
accuracy.  This  also  casts  doubt  on  the  assumed  species 
pdf  model’s  accuracy  for  computing  the  mean  species 
reaction  rate  (Equation  (1))  which  often  contains  higher 
moments  as  well.  This  deficiency  of  the  model  is  likely 
the  source  of  the  ineffectiveness  of  the  assumed  species 
pdf  formulation. 

VI.  CONCLUSIONS 

A  computational  study  was  undertaken  to  assess 
the  impact  of  turbulent-chemistry  interactions  on  the 
afterburning  and  afterburning  shutdown  characteristics  of  a 
generic  amine  booster.  Turbulent-chemistry  interactions 


were  accounted  for  using  an  assumed  probability  density 
function  based  method.  Analysis  of  the  simulation  results 
including  the  model  led  to  the  following  conclusions: 

1)  Turbulent-chemistry  interactions  modeled 
using  the  assumed  pdf  method  were  found  to  be  a  first 
order  effect  and  have  a  large  impact  on  plume  signatures  as 
afterburning  shutdown  was  approached.  The  pdf  model 
enhanced  plume  afterburning  and  delayed  the  onset  of 
shutdown  causing  the  plume  to  burn  to  higher  altitudes 
than  when  the  model  was  not  included.  However,  the  pdf 
model  did  not  change  the  character  of  the  shutdown  event. 
The  model  produced  a  gradual  type  of  shutdown  behavior 
characteristic  of  a  Damkohler  number  mechanism 
(Calhoon,  2000).  The  model  does  not  account  for  a  rapid 
shutdown  phenomenon  identified  by  Calhoon  (2000). 

2)  The  assumed  pdf  model  was  found  to  become 
increasingly  important  as  altitude  increased.  This  resulted 
from  the  increasing  trend  of  temperature  fluctuations 
within  the  shear  layer.  This  trend  was  found  to  be  a  result 
of  shock-shear  layer  interactions  within  the  plume. 

3)  The  calibration  constant  on  the  production  term 
in  the  energy  variance  equation  obtained  from  optimization 
using  low  speed  jet  data  was  found  to  be  inapplicable  to 
high  speed  plume  flows.  This  low  speed  calibration 
produced  excessively  high  temperature  fluctuations  which 
resulted  in  unrealistic  results.  Lower  values  of  the 
coefficient  produced  realistic  results  and  demonstrated  the 
importance  of  turbulent-chemistry  interaction  within  plume 
flows.  This  issue,  however,  highlights  the  need  for  a 
rigorous  validation  of  the  energy  variance  equation  for 
high  speed  flows  before  quantitative  predictions  can  be 
made. 

4)  The  species  fluctuation  aspect  of  the  assumed 
pdf  formulation  was  found  to  be  ineffective  in  the  present 
plume  flows.  Apparent  inaccuracies  in  the  evaluation  of 
the  chemical  source  term  in  the  scalar  energy  equation 
resulted  in  the  model  turning  itself  off.  This  same  trend 
has  been  observed  by  other  investigators  and  suggests  that 
a  reformulation  of  the  species  fluctuation  aspect  of  the 
model  is  warranted. 
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Fig.  1.  Illustration  of  the  behavior  of  plume  total  radiant 
intensity  as  a  function  of  altitude,  (a)  afterburning  regime, 
(b)  afterburning  shutdown  regime  (c)  post-afterburning 
shutdown  regime. 


T(K):  300  1250  Xcol:  0.0  0.1 


(a)  (b) 

Fig.  2.  Contours  of  (a)  temperature  and  (b)  C02  mole 
fraction  at  30  km  for  both  the  laminar  rate  (lam)  and 
temperature  pdf  (t  pdf)  models.  The  nozzle  exit  plane  is  at 
the  top  and  downstream  is  to  the  bottom.  The  transverse 
direction  has  been  scaled  by  a  factor  of  4  for  clarity. 


Fig.  3.  Comparison  of  laminar  and  temperature  pdf  rate 
model  predictions  of  centerline  mean  temperature  at  25, 
30  and  35  km  (a-  =  0  is  at  the  engine  nozzle  exit  plane). 
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Fig.  4.  Prediction  of  transverse  mean  CO:  mole  fraction  at 
x/xrefl  =  3  as  a  function  of  altitude  for  the  laminar  rate  and 
temperature  pdf  models. 


Fig.  5.  Prediction  of  transverse  mean  temperature  at  x/xreji 
=  3  as  a  function  of  altitude  for  the  laminar  rate  and 
temperature  pdf  models. 


Fig.  6.  Prediction  of  transverse  mean  OH  mole  fraction  at 
x/xreji  =  3  as  a  function  of  altitude  for  the  laminar  rate  and 
temperature  pdf  models. 


Fig.  7.  Temperature  intensity  predictions  from  the 
assumed  pdf  model  as  a  function  of  altitude  at  x/xreji  =  3. 


Fig.  8.  Transverse  variation  of  the  production,  dissipation 
and  dilatation  terms  in  Equ.  (12)  at  x/xreji  =  3  as  a  function 
of  altitude. 


Fig.  9.  Transverse  variation  of  the  sum  total  of  the 
production,  dissipation  and  dilatation  terms  in  Equ.  (12)  at 
x/xreii  =  3  as  a  function  of  altitude. 
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Fig.  10.  Source  station  radiation  variation  as  a  function  of 
altitude  at  zero  aspect  angle. 
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Fig.  11.  Comparison  of  laminar  rate  and  temperature  pdf 
model  predictions  of  source  total  radiant  intensity  as  a 
function  of  altitude  at  zero  aspect  angle. 
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Fig.  12.  Sensitivity  of  mean  temperature  to  the  calibration 
constant  Ci  in  Equ.  (12)  at  x/xre/i  =  3  for  30  km. 


Fig.  13.  Sensitivity  of  temperature  intensity  to  the 
calibration  constant  Ct  in  Equ.  (12)  at  x/xreji  =  3  for  30  km. 
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Fig.  14.  Prediction  of  transverse  mean  temperature  at 
x/xrefl  =  3  as  a  function  of  altitude  using  the  temperature 
pdf  only  (t  pdf)  and  the  temperature-species  pdf  (ty  pdf) 
formulations. 


Fig.  15.  Prediction  of  scalar  energy  from  Equ.  (13)  both 
with  and  without  the  chemical  source  term  included  at 

x/xrefl  =  3  and  30  km. 
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